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Abstract 

Within the framework of selfconsistent cranked Hartree-Fock- Bogoli- 
ubov theory(one-dimensional) we predict second backbend in the yrast 
line of 182 Os at I « 40ft, which is even sharper than the first one 
observed experimentally at / ~ 14ft. Around such a high spin the 
structure becomes multi-quasiparticle type, but the main source of 
this strong discontinuity is a sudden large alignment of z 13/ / 2 proton 
orbitals along the rotation axis followed soon by the alignment of j 15 / 2 
neutron orbitals. This leads to drastic structural changes at such high 
spins. When experimentally confirmed, this will be observed for the 
first time in this mass region, and will be at the highest spin so far. 



PACS numbers: 21.60.Jz, 21.10.Re, 21.60.Ev, 27.70.+q 
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During the latter half of the seventies and early eighties when the research 
activities in the high spin spectroscopy were rather at its peak, experimen- 
tally as well as theoretically, it was the one- dimensional cranking (rotation 
about a principal axis(x) perpendicular to the symmetry axis(z)) Hartree- 
Fock-Bogoliubov (CHFB) theory [0-01 which was most successful (besides 
its well-known shortcomings) to explain the underlying mechanism of the 
backbending features observed in the moment of inertia (or spin) versus ro- 
tational frequency plot of the yrast levels (level with the lowest energy for a 
given angular momentum). This is because of the fact that the mechanism of 
the alignment of single particle angular momenta along the rotation axis(the 
effect of collective rotation on single particle motion) is naturally present in 
this approach. 

After exhaustive activities in the rare-earth region, now the high spin 
structure work has spread to all the mass regions of the periodic table with 
the maximum spin value reached being of the order of I = 507i(it is 60h for the 
first superdeformed band in 152 Dy). Besides the yrast sequence, several side 
bands are observed in most of the nuclei. Recently very interesting features 
have been observed in some of the nuclei near the upper end of the rare-earth 



region, namely, in W — Os isotopes with mass number A ~ 180 || - [[13 
In some even - even isotopes side bands are found with high-K band heads 
(K being projection of I along the z-axis) very close to the yrast line with 
signature ( r = (— l) 7 ) symmetry broken, that is, even and odd spin states 
are connected by B(M1) transitions and the ratio B(M1)/ B(E2) is found 
to be large. Such states with mixed signature symmetry are interpreted 



as t bands |T3 arising due to the rotation of the nucleus about an axis 



tilted with respect to the principal axes of the quadrupole shaped deformed 
nucleus. Some of the states of these high-K bands decay to the yrast states, 
particularly in the band crossing region. So, it becomes natural to expect that 
the mutual interaction between these unperturbed bands should influence the 
backbending behaviour in the yrast sequence. Also in this mass region the 
Fermi surface lies in the high-m states of high-j orbitals(e.g. Ohu/2 for protons 
and 0zi3/ 2 for neutrons) which also gives some credence to such expectation 
that the yrast states, at least in the band crossing region, could be generated 
by tilted axis rotation. Thus, it has become an interesting and challenging 
problem to test if the usual explanation of backbending caused by the crossing 
of the s band ( low-m 213/2 neutron aligned band) still holds. This question 
is not yet resolved through a microscopic quantal many body calculation. 
It may be added that currently Onishi and his collaborators [15. 16 1 are 
attempting to perform a generator coordinate method (GCM) calculation for 
182 0s treating the tilting angles as generator coordinates. The good angular 
momentum projected GCM wave functions would be able to elucidate on the 
distribution of K components for a given value of the angular momentum. 

In some odd-A nuclei in this mass region it has recently been found that 
the high-K bands really cross the ground band producing backbend in the 
yrast sequence, e.g. l79 W |lTJ and 181 Re [12]. However, for even - even nuclei 



the general conclusion so far is that it is still the normal s band crossing that 
plays essential role to produce the backbend in these nuclei 0. Furthermore, 
in a recent microscopic theoreical analysis of the tilted axis rotation following 
the band mixing spin projected shell model approach |17| it is found that in 



the case of 178 > 180 ' 182 jy anc l 184 0s the backbend originates due to normal s 



band crossing; the tilted bands appear slightly above the yrast line. In view 
of these findings we thought of checking as to how the traditional CHFB ap- 
proach works for this mass region. Particularly for very high spins, / > 20h, 
the usual CHFB theory is expected to work well, and angular momentum 
projection is almost impracticable. Hence, we have carried out here a self- 
consistent CHFB calculation for 182 0s in an appropriate single paricle model 



space. For this nucleus presently the yrast line extends up to / = 34/2 || [13 
without a second backbend and the calculation has been performed for spins 
up to / = 50h. We find a strong backbend at around I = 40h. In the follow- 
ings we present our results and discussions on these. Also now on the spin 
values will be understood to be in the units of h. 

Since the CHFB theory 0, |4[] is well-known we will not give any de- 
tails here. For the Hamiltonian of the system we have used a pairing-plus- 
quadrupole model interaction. However, in addition, a hexadecapole term is 
also considered as these nuclei are expected to have large (negative) hexade- 



caploe deformation, /? 4 [13]. The total Hamilonian can be written as 



H = Ho ~ \ E XxEQ^(- i rQx-» -j E G T P}P T (1) 

A=2,4 A* r=p,n 

where x\ an d G T are the corresponding interaction strengths (in MeV), and 
the multipole moments (Q\fj) and the pairing operator (P) have standard 
forms 0, the radial part of the former being (r 2 /b 2 ) with b as the oscillator 
length parameter. Hq represents the spherical part of the hamiltonian with 
the single particle (s.p.) orbitals(assuming Z = 40, N = 70 core) 2si/2, 
14/2, O07/2, 009/2, Oiis/a, I/7/2, 0/i 9/2 , and 0h u/2 for protons, and 2p 1/2 , 2p 3/2 , 
I/5/2, I/7/2, O/ig/2, O/in/2, Oj'15/2, 1(79/2, 0iu/ 2 , and 0ii3 /2 for neutrons. The 



s.p. energies are as given in table. 1 of Ref. [15j except that of Oj'15/2 which is 
7.179 MeV. These energies are essentially the Nilsson spherical s.p. energies 
for this mass region. 

As indicated above, Eq. ([]]) actually represents two Hamiltonians as 
far numerical calculations are concerned here. One is as it is, and in the 
other the hexadecapole term is dropped. Correspondingly we have two sets 
of ineraction strengths which can reproduce more or less the ground state 
intrinsic shape parameters of 182 0s. Finally taking principal x-axis as the 
cranking axis the CHFB eigenvalue equations are solved by diagonalization 
with the usual particle number and angular momentum constraints [fj. |J]. 

We have also computed rotational g Factors using the standard cranking 
expression || 

9I = fa) /y/ 1(1 + 1) (2) 
where fi x is the x-component of the magnetic moment operator 

/* = 0i£j*(i) + (0« ( 3 ) 

i i 

with gi — 1 and g s = 5.586 for protons and g% = and g s = —3.826 for 
neutrons. In numerical computation the values of g s are attenuated by a 
factor of 0.6 g, . 



The CHFB equations are solved selfconsistently in terms of seven col- 
lective variables when the hexadecapole(Q4) term is also present in the 
Hamiltonian. These are the pairing gaps, A p and A n , for protons and neu- 
trons, respectively and deformation parameters = q\-^ = {Q\fj,) with 
H — 0, 2 for A = 2 and fi = 0, 2, 4 for A = 4. The usual deformation param- 
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eters ( (3, 7 ) and [3± are defined through the relations, hu>oP cos 7 = 
frcuo sin 7/ = X2?22 and = X4<Z40> where ftco> = 4.1/ A 1 / 3 MeV. 

When only quadrupole interaction is considered then the degree of free- 
dom is reduced to four. As already mentioned, the interaction strengths 
are chosen such that the values of the ground state shape parameters (see 
Table I) are approximately reproduced || [18[]. The values of the pairing 



gap parameters are decided by looking at the experimental odd-even mass 
differences. It may be emphasized that after fixing the interaction strength 
parameters at this stage there are no free parameters in the theory. In the 
followings, for the sake of brevity, when the quqdrupole + hexadecapole both 
the interaction terms are considered , then the corresponding results will be 
indicated by the symbol "Q4". But if only quadrupole term is considered, 
then it will accordingly be indicated by "Q2". 

Now we can discuss some of our main results. In Table I we have listed 
the values of the shape parameters at a few angular momentum values, the 
dependence on / being most striking around / = 40. In both the cases, with 
and without the Q4 term, A n goes to zero at / = 14. On the other hand A p 
vanishes at / = 30 in presence of the Q4 term, and at / = 26 without it which 
implies a somewhat stronger proton pairing correlation in the former case for 
the spin range of about I = 16 — 28. We find that when only quadrupole 
force is considered the proton pairing recovers for / = 34 (A p = 0.214 MeV) 
to / = 40 (A p = 0.239 MeV). In the other case it recovers only at one spin, 
/ = 34 with a small value of A p = 0.174 MeV. In both the cases there is a 
correlation between 7 acquiring a negative value and an increase in the value 
of /3. Around spin / = 40 7 changes sign from positive to negative by a quite 
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sizeable amount, and associated with it the value of (3 increases by about 
20% (stretching). At the same time the value of /3 4 shows a sudden decrease. 
However, it may be pointed out that now all the three components of q^, 
/i = 0,2,4 become of similar magnitude: for instance, at / = 38 these are 
-8.05, -1.09 and 1.08 which become -4.10, -2.36 and 2.56 (all in units of b 2 ) 
at / = 40 for /i = 0, 2 and 4, respectively. 

In Fig. [I] we display a backbending(BB) plot of spin versus rotational fre- 
quency (u) where for the experimental case u>i = \(Ei — Ej^ 2 )- As indicated 
on the top right corner of the figure, the three curves correspond to the exper- 
imental data, with hexadecapole(Q4) and without hexadecapole(Q2) terms 
in the Hamiltonian of Eq. (H). We notice that the first backbend is actu- 
ally not well reproduced, though an upbend is produced at more or less the 
correct frequency, and the inclusion of the hexadecapole degrees of freedom 
helps in the right direction. However, the alignment of O213/2 neutron orbitals 
is quite pronounced as can be seen in the next Fig. 0, where contributions 
from a few important orbitals to the total angular momentum are shown 
for the "Q4" case. At / = 14 the contribution of m 13 / 2 orbitals (mainly 
m = 7/2,9/2 components) is about 9h, close to the experimental estimate 
of 10h [|J. However, the alignment is not sudden enough around this spin to 
cause a sharp BB. Hence, we may conclude that some extra mechanism is, 
perhaps, needed to obtain a sharp first BB. We are trying to perform angular 
momentum projection, including K-mixing, on CHFB wave functions. If this 
also fails, then, perhaps, tilting mechanism is the only explanation. 

However, the main interesting result here is the appearance of a second 
BB near / = 40 as seen in Fig. |I[ The small backward kink in the " Q4" curve 
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at I = 36 is a genuine one, that is, it is not due to some numerical inaccuracy 
etc. The sharp BB is, of course, due to a sudden large and coherent alignment 
of the pi 13/2, phg/2 and nj'15/2 orbitals (mainly m = 1/2 and 3/2 components) 
as clearly seen in Fig. |2|. At / = 40 the contribution of ra 13 / 2 suddenly drops 
by about 6 units in one step. The nhg/2 orbitals contribute at all spins, in 
almost a gradual manner, whereas p/111/2 orbitals start contributing at very 
high spins through high-m components. Thus, the structure near / = 40 is 
very interesting. There is a quite sizeable stretching of /3, and 7 acquires a 
negative value. That is in this BB region collectivity increases, rather than 
showing a decrease, as is often observed in the first band crossing region 
[ |j~9| , p0| . For the spin region / = 40 — 50 the intrinsic structure remains 
essentially unchanged, which may be seen as a second minimum in a shell 
correction calculation, and levels for / > 40 may be interpreted as rotational 
levels in the second well. The increase in (3 can be understood as due to an 
enhanced magnitude of the quadrupole matrix elements of the low-m, high-j 
(aligned) neutron and proton orbitals near the Fermi surface. 

In Fig. ^ is displayed a variation of the g Factors with spin. The actual 
value of g(I = 2) is 0.245 for the "Q4" case and 0.230 for the "Q2" case 
which appear to be reasonable in view of gn ■ = 0.27 used in Ref. ||. In both 
the cases the ratio g(I)/g(2) drops sharply to a minimum at I = 14, though 
in " Q4" case the real minimum is at I = 20 with a slightly lower value. The 
sharp rise at I = 40 is a clear indication of the large alignment of the proton 
orbitals. Beyond this the alignment of nj'15/2 orbitals stops the further rise. 
In the intermediate spin region (/ = 20 — 30) relatively larger magnitude of 
the alignment of the large-m components of phn/2 make the g Factors higher 
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for the "Q2" case. 

Finally we would like to make some additional remarks. Since at / = 40 
the pairing has collapsed in our calculation, one may think that in a particle 
number projected treatment the position of BB may get shifted or become 
much less dramatic. In order to check for this a calculation was performed 
for / > 30 with a fixed value of A p and A n at about half of their values in 
the ground state. Then it is found that 7 changes sign between / = 32 to 34 
with all the characteristics as noted above. Thus, the second BB seen here 
seems very much genuine. We also notice a small favourable trend looking 
at the difference of the experimental 7 - ray energies: E y (32) - E y (30) = 78 
keV, and £ 7 (34) - £ 7 (32) = 65 keV. 

In conclusion, through a selfconsistent CHFB calculation, which is very 
reliable at high spins (/ > 20), we have obtained a clear case of sharp second 
backbend in 182 0s near / = 40. This is caused by a large coherent alignment 
of low-m piu/2, phg/2 and nj'15/2 orbitals. In this spin region there is a sub- 
stantial change of structure within a couple of units of angular momentum. 
For instance, 7 goes positive to negative, with a change of 12° or more in one 
step, with an associated increase of f3 by about 20%(stretching). However, 
the first backbend is not well reproduced, and, as discussed above, we are 
working on it. We have also studied the variation of g Factors as a function 
of spin which essentially shows, in a gross manner, the alignment pattern of 
neutron and proton single particle orbitals. 

Particularly for 182 0s the levels are already known up to / = 34, so we 
hope that our prediction will produce enough excitement in experimentalists 
to put efforts to study the interesting features in the spin / = 40 region. 
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Table 1: Intrinsic shape parameters as a function of spin, I (actually 
\Jl{I + 1))- For the ground state 7 = and (3,04, A p and A n are, respec- 
tively 0.228, -0.038, 0.871 MeV, and 0.879 MeV in case "Q4", and 0.229, 0.0, 
0.872 MeV and 0.886 MeV in case "Q2". At I = 10 A n = 0.436 MeV in the 
former case and 0.432 MeV in the latter case. It goes to zero for I > 14. 
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Figure Captions: 



Figure 1: Backbend plot for 182 0s showing the variation of angular momen- 
tum with rotational frequency, uo. As indicated on the figure the curve with 
solid line corresponds to the experimental data with a>j = \{Ei — Ei-%). The 
long dashed curve, labeled "Q4" corresponds to the case when quadrupole 
and hexadecapole both the interaction terms are present in the Hamiltonian 
flip. The short dashed curve, labeled "Q2" indicates that only quadrupole 
interaction is considered. 

Figure 2: Alignment plot. Contributions of a few important orbitals to 
the total spin is displayed denoting it as a single particle (S. P.) contribution. 
These are corresponding to the "Q4" case only. The type of orbitals is indi- 
cated for each curve, where p and n indicate proton and neutron, respectively. 

Figure 3: Variation of the g Factor ratio g(I)/g(2) as a function of spin 
for both the cases, that is, with the inclusion of hexadecapole term(Q4), and 
without it(Q2). 
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Alignment plot for Os-182 
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